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Abstract 

We explore the momentum and velocity dependent elastic scattering between the dark matter (DM) par- 
ticles and the nuclei in detectors and the Sun. In terms of the non-relativistic effective theory, we phe- 
nomenologically discuss ten kinds of momentum and velocity dependent DM-nucleus interactions and re- 
calculate the corresponding upper limits on the spin-independent DM-nucleon scattering cross section from 
the current direct detection experiments. The DM solar capture rate is calculated for each interaction. Our 
numerical results show that the momentum and velocity dependent cases can give larger solar capture rate 
than the usual contact interaction case for almost the whole parameter space. On the other hand, we deduce 
the Super-Kamiokande's constraints on the solar capture rate for eight typical DM annihilation channels. 
In contrast to the usual contact interaction, the Super-Kamiokande and IceCube experiments can give more 
stringent limits on the DM-nucleon elastic scattering cross section than the current direct detection experi- 
ments for several momentum and velocity dependent DM-nucleus interactions. In addition, we investigate 
the mediator mass's effect on the DM elastic scattering cross section and solar capture rate. 
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I. INTRODUCTION 



The existence of dark matter (DM) is by now well confirmed |2j. The recent cosmological 
observations have helped to establish the concordance cosmological model where the present Uni- 
verse consists of about 68.3% dark energy, 26.8% dark matter and 4.9% atoms yfl. Understanding 
the nature of dark matter is one of the most challenging problems in particle physics and cosmol- 
ogy. The DM direct detection experiments may observe the elastic scattering of DM particles with 
nuclei in detectors. Current and future DM direct search experiments may constrain or discover 
the DM for its mass m D and elastic scattering cross section cr„ with nucleon. As well as in the DM 
direct detection, the DM particles can also elastically scatter with nuclei in the Sun. Then they 
may lose most of their energy and are trapped by the Sun [1]. It is clear that the DM solar capture 
rate C© is related to the DM-nucleon elastic scattering cross section <x„. Due to the interactions 
of the DM annihilation products in the Sun, only the neutrino can escape from the Sun and reach 
the Earth. Therefore, the water Cherenkov detector Super- Kamiokande (SK) [4] and the neutrino 
telescope IceCube (IC) oils] can also give the information about m D and <x„ through detecting the 
neutrino induced upgoing muons. 

The current experimental results about cr n are based on the standard DM-nucleus contact inter- 
action which is independent of the transferred momentum q and the DM velocity v. In fact, many 
DM scenarios can induce the momentum and velocity dependent DM-nucleus interactions. For 
example, the differential scattering cross section of a long-range interaction will contain a factor 
(q 2 + rn^y 2 with being the mass of a light mediator QVL ISQ. It is worthwhile to stress that 
the current experimental results about <x„ must be recalculated for the momentum and velocity 
dependent DM-nucleus interactions. In view of this feature, many authors have recently used the 
momentum and velocity dependent DM-nucleus interactions to reconcile or improve the tension 
between the DAMA annual modulation signal and other null observations Il8l-ll2ll. The new upper 
limit on <j n can directly affect the maximal C Q . On the other hand, we have to recalculate C© for 
a fixed cr n when the DM-nucleus interaction is dependent on the momentum and velocity. For the 
usual contact interaction, the current direct search experiment XENON100 [13;] provides a more 
stringent limit on spin-independent (SI) <x„ than the Super-Kamiokande and IceCube experiments 
when m D > 10 GeV fl-Q]. We do not know whether this conclusion still holds for the momen- 
tum and velocity dependent DM-nucleus interactions. It is very necessary for us to systematically 
investigate the momentum and velocity dependent DM elastic scattering in detectors and the Sun. 
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In this paper, we will explore the momentum and velocity dependent DM-nucleus interactions 
and discuss their effects on the SI cr„ and the DM solar capture rate C Q . New upper limits on 



cr„ from the XENON100 EM and XENON10 EM, and the corresponding maximal C Q will be 



calculated for these interactions. On the other hand, we shall deduce the constraints on C Q from 
the latest Super- Kamiokande results for eight typical DM annihilation channels. In addition, the 
mediator mass's effect on <j n and C Q will also be analyzed. This paper is organized as follows: 
In Sec. II, we outline the main features of the momentum and velocity dependent DM-nucleus 
interactions in direct detection experiments, and derive the corresponding upper limits on cr n . In 
Sec. Ill, we numerically calculate C Q for these interactions and give the general constraints on C Q 
from the Super-Kamiokande and IceCube. In Sec. IV, we discuss the mediator mass's effect on cr n 
and C Q . Finally, some discussions and conclusions are given in Sec. V. 

H. DARK MATTER DIRECT DETECTION 
A. DM event rate 



The event rate R of a DM detector in the direct detection experiments can be written as 

r = Nt p™ r d^ dER p v/(v)J 3 V5 

J dt R Jv 



= N, 



Pdm nA 2 m N o- n 



/pi p' 'max 

F 2 N (q)dE R I dcosO I vf(v)F% M (q, v)dv , (1) 



where N T is the number of target nucleus in the detector, p DM = 0.3GeVcm 3 is the local DM 
density, mo is the DM mass. For the DM-nucleus differential scattering cross section dcr N /dE R , 



we have taken the following form II 1 0(1 



dcr N A2 m N o- n 2 2 

drr A ^ FNiq)F ™ (q ' v) ' (2) 

where A is the mass number of target nucleus, cr n is the DM-nucleon scattering cross section. 
In Eq. ©, we have required that the proton and neutron have the same contribution. The DM- 
nucleon reduced mass is given by jx n = m D m n /(m D + m n ) where m n is the nucleon mass. The 
recoil energy E R is related to the transferred momentum q and the target nucleus mass m N through 
q 2 = 2m N E R . The DM velocity distribution function f(y) in the galactic frame is usually assumed 
to be the Maxwell-Boltzmann distribution with velocity dispersion vq = 220 km/sec, truncated at 
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the galactic escape velocity v esc = 544 km/sec. In the Earth's rest frame, we can derive 

W = r4^~ (l7+l7f)2/v ° . ( 3 ) 

(nv 2 ) 3 ' 2 

where v is the DM velocity with respect to the Earth and v e « v© = vo + 12km/sec is the Earth's 
speed relative to the galactic halo. It is worthwhile to stress that the contribution of the Earth's 
orbit velocity to v e has been neglected since we do not focus on the annual modulation. With the 
help of |v + v e \ < v esC5 one can obtain the maximum DM velocity 

Vmax = tJvI sc - V 2 e + v\ COS 6 2 - V e COS Q , (4) 

where 6 is the angle between v and v e . For a given recoil energy E R , one can easily derive the 
minimum DM velocity 

^2m N E R 

Vmin = » , (5) 

2/ijv 

where ///v = m D m N /(m D +m N ) is the DM-nucleus reduced mass. For the nuclear form factor Fj,(q), 



we use the Helm form factor 



3j'i(^i) 
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with 7?! = ^c 2 + \n 2 a 2 -5s 2 and c - 1.23A 1/3 - 0.60 fm [16]. Here we take s - 0.9 fm and 
a - 0.52 fm IU6I1 . j\(x) = sinx/x 2 - cosx/x is a spherical Bessel function of the first kind. For 



the usual contact interaction, the DM form factor F^ M (q, v) = 1 is independent of the transferred 
momentum q and the DM relative velocity v. In this paper, we shall focus on some momentum 
and velocity dependent DM form factors and discuss their effects on the DM direct detection cross 
section and the DM solar capture rate. 



B. Momentum and velocity dependent DM form factors 

Usually, one can build some DM models and exactly calculate the DM direct detection cross 
section. On the other hand, the DM-nucleus interaction can be generally constructed from 16 



model-independent operators in the non-relativistic (NR) limit I117U 18I1. Any other scalar operators 
involving at least one of the two spins can be expressed as a linear combination of the 16 indepen- 
dent operators with SI coefficients that may depend on q 2 and V* 2 = (v-q/(2[i N )) 2 = v 2 -q 2 /(4fj^ / ). 
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TABLE I: The momentum and velocity dependent DM form factors for the heavy and light mediator mass 



scenarios with V 2 = v 2 - q z /(4fii,), q re f = 100 MeV and V re f = vq. 



It is convenient for us to phenomenologically analyze the momentum and velocity dependent DM- 
nucleus interactions from these N R op erators. Here we only focus on the following four SI NR 



operators in the momentum space [117 



18b: 



2 
3 
4 



1 , 

is D ■ q , 
s D -V , 
i? D ■ (V 4 x 



(7) 



Considering the possible contributions of q 2 or V 2 in the coefficients, we phenomenologically 
discuss five kinds of momentum and velocity dependent DM form factors F 2 M (g, v) up to q and 
V quartic terms in the amplitude squared \M\ 2 . The five DM form factors and the usual contact 
interaction case have been listed in the third column of Table H Since the transferred momentum q 
in many direct detection experiments is order of 100 MeV, we take q K f = 100 MeV as the reference 
transferred momentum to normalize q. Similarly, we use V K f = v Q to normalize V = ^jv 2 - . 
Here we have assumed that the mass of mediator between DM particles and quarks is far larger 
than the transferred momentum q, namely :» q 2 . If <sc q 2 , Fp M (^, v) should contain the 
factor 1 lq 4 which comes from the squared propagator (q 2 + m 2 ^)' 2 . For the light mediator mass 
scenario, the corresponding 6 kinds of F^ M (q, v) cases have been listed in the sixth column of 
TablelH In Sec. IV, we shall discuss the m 2 ~ 0(q 2 ) scenario through varying m^. 
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C. New upper limits on cr n 



In this paper, we do not try to reconcile the tension between the DAMA annual modulation 
signal and other direct detection exclusions by use of the momentum and velocity dependent DM- 
nucleus interactions. Here we only focus on the null observations and the corresponding upper 
limits on <r n which are relevant to the maximal DM solar capture rate. Currently, the most stringent 



limit on a n comes from XENON 100 lij and XENON 10 yj]. It should be mentioned that this 



limit is only valid for the usual contact interaction, namely F^ M (q, v) = 1 . For the momentum 
and velocity dependent F^ M (q, v), we should recalculate their limits from the reported results of 
XENON 100 and XENON 10. 

The recoil energy window of the DM search region in the XENON 100 is chosen between 
3 ~ 30 photoelectrons (PE), corresponding to 6.6 keV < E R < 43.3 keV. The relation of E R and 



PE number S 1 is given by IU3I1 



S \(E R ) = 3.73 PE x E R x £ eS , (8) 

where X e ff is the scintillation efficiency which has been measured above 3 keV. The X e ff 
parametrization can be found in Ref. Ill9ll . Here we assume that the produced PE number of a 
nucleus recoil event satisfies the poissonian distribution and Eq. © denotes the mean value. In 
this case, the event with E R < 6.6 keV will have a non vanishing probability to generate a S 1 
signal above 3 PE. For the new lower threshold of E R , we take E R > 3.0 keV which can pass the 



ionization yield S 2 cut 11 



The search recoil energy range of XENON 10 is 1 .4 keV < E R < 10.0 keV [UJ. For 4 GeV < 
m D < 20 GeV, one can always find some parameter space among 1.4 keV < E R < 10.0 keV to 
satisfy v m i n < v max . Therefore, we directly input 1.4 keV < E R < 10.0 keV into Eq. (QQ) for the 
XENON10 analysis. Note that the upper limit with vo = 230 km/sec and v esc = 600 km/sec in Ref. 



has been replaced by the corresponding limit with vo = 220 km/sec and v esc = 544 km/sec in 



the following parts. 

Requiring the same event rate R for different F^ M (q,v), we deduce new bounds about cr n for 
each Fl M (q, v) from the F^ M (q, v) = 1 case (the reported limits of XENON100 and XENON10). 
Our numerical results have been shown in the left panel of Fig. \T\ The same color solid and dashed 
lines describe the heavy and the corresponding light mediator mass scenarios, respectively. In Fig. 
CD the number 1 denotes the F^ M (q, v) = 1 case, V 2 q~ 2 denotes the F^ M {q, v) = V 2 q 2 ef /(V^q 2 ) 
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FIG. 1: The new upper limits on cr n for different F^ M (q, v) with g re f = 100 MeV; V re f = vo (left panel) and 
q Te{ = 10 VlO MeV; V ief = 2v (right panel) from the XENON100 and XENON10. 



case, and so on. It is meaningless for us to compare different lines since these limits are dependent 
on q ie { and V re f. For illustration, we plot our numerical results with q re{ = 10 VlO MeV and 
V re f = 2vo in the right panel of Fig. \T\ Some kinks around m D = 8 GeV arise from the different 
slopes of the predicted limits of XENON 100 and XENON 10. It is should be mentioned that the 
new upper bound from the XENON 100 and XENON 10 is still the most stringent limit for each 



F^ M (q, v) when we recalculate other experimental results [12111 . 



III. DARK MATTER SOLAR CAPTURE 



When the halo DM particles elastically scatter with nuclei in the Sun, they may lose most 
of their energy and are trapped by the Sun [1]. On the other hand, the DM annihilation in the 
Sun depletes the DM population. The evolution of the DM number N in the Sun is given by the 
following equation [12211 : 



N = C Q - C a N 2 , 



(9) 



where the dot denotes differentiation with respect to time. The DM solar capture rate C Q in Eq. © 
is proportional to the DM-nucleon scattering cross section cr n . In the next subsection, we shall give 
the exact formulas to calculate C . The last term C A N 2 in Eq. © controls the DM annihilation rate 
in the Sun. The coefficient C A depends on the thermal-average of the annihilation cross section 
times the relative velocity (<xv) and the DM distribution in the Sun. To a good approximation, one 
can obtain C A = (<xv)/V e ff, where V eff = 5.8 x 10 30 cm 3 (lGeV/m £) ) 3/2 is the effective volume of 



the core of the Sun Q22 



2311 . In Eg. ©, we have neglected the evaporation effect since this effect 



is very small when m D > 4 Ge V 11241. 12511. One can easily solve the evolution equation and derive 
the DM solar annihilation rate Il22ll 



T A = ic tanh 2 a o VC G C A ) , 



(10) 



where t e - 4.5 Gyr is the age of the solar system. If t Q VQdCa :» 1, the DM annihilation rate 
reaches equilibrium with the DM capture rate. In this case, we derive the maximal DM annihilation 
rate Y A = C Q /2. It is clear that the DM annihilation signals from the Sun are entirely determined 
by C Q . 



A. DM solar capture rate and annihilation rate 

By use of the DM angular momentum conservation in the solar gravitational field, one can 



obtain the following DM capture rate C Q 112311 : 

C Q = J]J 4nr 2 dr J ywn N »J 3 « (11) 

with 

fto = r4^~ (lt+lt0)2/v ° . a 2 ) 

(nv 2 y /2 

where f(u) is the DM velocity distribution, u is the DM velocity at infinity with respect to the 
Sun's rest frame, v Q = v + 12 km/sec is the Sun's speed relative to the galactic halo. 0. N .{cS) is the 
rate per unit time at which a DM particle with the incident velocity co scatters to an orbit within 
the Jupiter's orbit. 0. Ni (cS) is given by 

Cl N , (co) = n N . (r)cr Ni (co)cop DM /m D , (13) 

where n^^r) and oj(r) = -yju 2 + Vg SC (r) are the number density of element and the DM incident 
velocity at radius r inside the Sun, respectively. The escape velocity v esc (r)from the Sun at the 
radius r can be approximately written as Vg SC (r) = v 2 c - {v 2 c - v 2 )M(r)/M Q I126I1 . where v c = 1354 
km/sec and v v = 795 km/sec are the escape velocity at the Sun's center and surface, respectively. 
M = 1.989 x 10 33 g is the solar mass and M{r) is the mass within the radius r. <r Ni {pS) in Eq. (TOT) 
is the scattering cross section between a stationary target nucleus N, in the Sun and an incident 
DM particle with velocity co. The non-relativistic effective theory allows us to express cr^a/) as 

A 2 cr n r 2 ^ 



cr Ni (to) 



5 I F 2 Ni (q)Fl M (q,co)qdq , (14) 
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where q m m = -\Jm D m Nj [u 2 + Vg SC (r = 5.2AU)] is the minimum transferred momentum needed for 
capture and v e ^Xr = 5.2AU) = 18.5 km/s denotes the DM escape velocity from the Sun at the 



Jupiter's orbit [17, 



220. 
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FIG. 2: The predicted maximal DM solar capture rates C Q for the heavy (left panel) and light (right panel) 
mediator mass scenarios. 



In our calculation, we sum over the following elements in the Sun: l H, 4 He, 12 C, 14 N, 16 0, 17 0, 
Ne, Mg, Si, S and Fe. The number densities n N .(r) of these elements and M(r) can be obtained 
from the calculation of the standard solar model (SSM). Here we employ the SSM GS98 fl28h to 
calculate the DM solar capture rate C G in Eq. (fTTj) with the help of <r„ in the left panel of Fig. 
[D Our numerical results have been shown in Fig. |2] We find that for almost whole of the m D 
parameter space the predicted C Q from the standard contact interaction is smaller than those from 
the momentum and velocity dependent DM form factor cases. This means that the momentum and 
velocity dependent DM form factor cases can give larger DM annihilation signals than that from 
the usual contact interaction case. The same color solid and dashed lines in Fig. [2] describe the 
heavy and the corresponding light mediator mass scenarios, respectively. The light mediator mass 
scenario can usually produce the larger C Q than the corresponding heavy mediator mass scenario. 
However one can derive the opposite conclusion for the q 4 and 1 cases. It should be mentioned 
that our numerical results in Fig. [2] are independent of q ref and V iei . 
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B. Constraints from the Super-Kamiokande and IceCube 



Due to the interactions of the DM annihilation products in the Sun, only the neutrino can escape 
from the Sun and reach the Earth. For the given DM mass and DM annihilation channel a, the 
differential muon neutrino flux at the surface of the Earth from per DM pair annihilation in the 
Sun can be written as 



dE V)i AnR\ s dE v 



(15) 



where R ES = 1.496 x 10 13 cm is the Earth-Sun distance. dN" JdE Vfi denotes the energy distribution 
of neutrinos at the surface of the Earth produced by the final state a through hadronization and 
decay processes in the core of the Sun. It should be mentioned that some produced particles, such 
as B mesons and muons, can lose a part of energy or the total energy before they decay due to 
their interactions in the Sun. In addition, we should consider the neutrino interactions in the Sun 



and neutrino oscillations. In this paper, we use the program package WimpSim [29] to calculate 



dN"JdE Vfi with the following neutrino oscillation parameters Q30 



SID: 



sin 2 n = 0.32, sin 2 # 23 = 0.49, sin 2 13 = 0.026, 6 = 0.83tt, 

Am 2 2l = 7.62 x lO^eV 2 , Am 2 3i = 2.53 x lO^eV 2 . (16) 

In addition, we should also calculate the differential muon anti-neutrino flux which can be evalu- 
ated by an equation similar to Eq. (fT5l) . 

These high energy neutrinos interact with the Earth rock or ice to produce upgoing muons 
which may be detected by the water Cherenkov detector Super-Kamiokande [4] and the neutrino 
telescope IceCube [|5|,|6|]. Due to the produced muons scattered from the primary neutrino direction 
and the multiple Coulomb scattering of muons on route to the detector, the final directions of 
muons are spread. For 10 GeV < m D < 1000 GeV, the cone half-angle which contains more 
than 90% of the expected event numbers ranges from 6° to 30° for the Super-Kamiokande when 
we assume the bb annihilation channel. The cone half-angles will be smaller for the other DM 
annihilation channels considered in this paper with the same DM mass. In terms of the results of 
cone half-angle 6 in Tables 1 and 2 of Ref. [|4J], we conservatively take some reasonable 6 for other 
DM annihilation channels and several representive DM masses as shown in Table UH 

The neutrino induced upgoing muon events in the Super-Kamiokande can be divided into three 
categories: stoping, non-showering through-going and showering through-going [4]. The fraction 
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TABLE II: The relevant parameter summary to calculate the Super-Kamiokande constraints on Ya for dif- 
ferent DM annihilation channels and masses. The units of mn and <pu are GeV and 10 _15 cm _2 sec _1 . 



of each upgoing muon category as a function of parent neutrino energy E v has been shown in Fig. 
2 of Ref. |4|]. Then we use dN"JdE Vii to calculate the fraction of each category F' as listed in 
Table HB Once F l is obtained, the 90% confidence level (CL) upper Poissonian limit N90 can be 
derived through the following formulas D4|] : 



J -CO 



90% 



and 



/=i 



(17) 



(18) 
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where v s is the expected real signal. The number of observed events of each category n l ohs and the 
expected background of each category n' BG for different DM masses and cone half-angles can be 
found in Tables 1 and 2 of Ref. Jj. With the help of Eqs. (QJ]) and CDS), we estimate the 90% CL 
upper Poissonian limit on the number of upgoing muon events N 90 and the corresponding 90% CL 
upper Poissonian limit of upgoing muon flux 0^ = N 90 /(l-67 x 10 15 cm 2 sec) as shown in TableUB 

->27 
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CD 
CO 

CN 

O 
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FIG. 3: The current Super-Kamiokande and IceCube constraints on C©/2 with the assumption Ta = C©/2 
and the predicted maximal DM solar capture rates C Q /2 from Fig. [2]for different F 2 M (g, v). The black solid 
line describes the equilibrium condition for (o~v) « 3.0 x 10~ 26 cm 3 sec -1 . 



With the help of Eq. (26) in Ref. I32Q, we numerically calculate the neutrino induced muon flux 
from per DM pair annihilation in the Sun. Then we directly derive the Super-Kamiokande con- 
straints on Y A from the 0^ values as listed in TableUIl In Fig. [31 we plot these results with the dotted 
lines and the predicted maximal DM solar capture rates C©/2 from Fig. [2] for different F^ M (q, v). 
It should be mentioned that Y A = C Q /2 has been assumed in Fig. [3] As shown in Eq. (fTOb . the 
assumption Y A = C Q /2 holds if t e VOjCa !• F° r me usual 5-wave thermally averaged annihi- 
lation cross section (crv) « 3.0 x 10 _26 cm 3 sec 1 deduced from the DM relic density, we find that 
t e VCqCI > 3.0 (namely tanh 2 [fe ^JC^C~a\ > 0.99) requires C Q /2 > 4.3 x 10 22 /(m D /lGeV) 3/2 sec- 1 
which has been plotted in Fig. [3] with the black solid line. Therefore the predicted C Q /2 above 
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this line in Fig. [3] will favor the assumption Y A = C Q /2. In addition to the Super-Kamiokande 
experiment, the IceCube collaboration has also reported the upper limits on the DM annihilation 
rate Y A for the bb and W + W~ (r + -r below m D = 80.4 GeV) channels in Table I of Ref. |5|]. We plot 
these results with the dash-dotted lines in Fig. |3j The IC&6(W + W~)* line shows the expected 180 
days sensitivity of the completed IceCube detector flg]. It is shown that the upper limits on C Q (cr„) 
from the Super-Kamiokande and IceCube are weaker than those from the current direct detection 
experiments for the usual SI DM-nucleus interaction. However, our numerical results in Fig. |3] 
clearly show the Super-Kamiokande and IceCube may give more stringent constraints than the 
XENON 100 experiment for several momentum and velocity dependent DM-nucleus interactions 
with m D > 10 GeV. 

In Fig. [3J one may find that the Super-Kamiokande experiment can significantly constrain the 
low mass DM for all DM form factors in Table [I] when the DM particles dominantly annihilate into 
neutrino pairs or t + t . If m D > 20 GeV, both Super-Kamiokande and IceCube can not constrain 
any momentum and velocity dependent case except for the V 4 q~ 4 case, when the annihilation 
channel is the bb, and the l,q 2 ,q~ 2 ,q 4 , V 2 , V 2 q 2 and V 2 q~ 2 cases for any annihilation channel. 
Except for the bb channel, other annihilation channels can be used to significantly constrain the 
V 4 q~ 4 and V 2 q~ 4 cases. For the W + W~ channel, the IceCube gives the stronger constraint than the 
Super-Kamiokande when m D > 100 GeV. The future IceCube result IC&6(W + W~)* has ability to 
constrain the q~ 4 , V 2 , V 4 , V 2 q 2 and V 2 q~ 2 cases with m D > 200 GeV. Since C Q is proportional to cr n , 
the upper limits on C Q in Fig. [3] will move downward if <x„ in Fig. [Qbecomes smaller. The Super- 
Kamiokande experiment can still constrain the momentum and velocity dependent DM-nucleus 
interactions for the low DM mass region even if cr„ is reduced by 2 orders. 



IV. THE MEDIATOR MASS EFFECT ON tr„ AND C Q 

In Sec. Ill BL we have taken two extreme scenarios for the mediator mass mf m 2 . :» q 2 and 
m 2 ^ <s q 2 . Here we shall consider the ~ 0(q 2 ) scenario and discuss the m$ effect on cr n and 
C Q . In this case, the momentum and velocity dependent DM form factors are relevant to m^. It is 
found that the dependent DM form factors F^ M (q, v, m^) can be written by the product of the 
third column of Table Hand a factor (q 2 e{ + m 2 ^) 2 /(q 2 + m^) 2 . The two DM form factors in each 
row of Table H are two extreme cases of the dependent DM form factor F^ M (^, v, m^). For 
example, one can easily obtain F 2 M (q, v) = q 2 /q 2 et - with m 2 » q 2 , q 2 and F^q, v) = qL/q 2 with 
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m 2 «; q 2 ,q 2 ei from F 2 M (g, v, m ) = (q 2 / q 2 ef )(q 2 ei + mj) 2 /(q 2 + m 2 ) 2 . Therefore, we have 6 kinds 
of m<£ dependent DM form factors F^ M (q, v, m$). Here we use lm^, q 2 m,p, V 2 m,p, q A m^, V 4 m^ and 
V 2 q 2 m l p to express them. 




10' 10" 1 
m (GeV) 



10 10 

m (GeV) 



FIG. 4: The predicted <x„ and C Q as a function of for m D = 10 GeV and m D = 100 GeV. The parameters 
Ra- n and R Ce are defined as R an = cr n /o- n (mj «; q 2 ,q 2 mi ) md R Ce - C Q /C Q (m 2 «; q 2 ,q 2 ei ). The solid and 
dashed lines describe mo = 10 GeV and mp - 100 GeV cases, respectively. 



Using the above 6 dependent DM form factors, we calculate cr n and C for two represen- 
tative DM masses: m D = 10 GeV and m D = 100 GeV. Our numerical results have been shown 
in Fig. |H The parameters R an and R Cq in Fig. |4]are defined as = o-Jo- n (m\ «; q 2 ,q 2 e{ ) and 
Rc B = C /C (m 2 q 2 , q 2 ef ). a n {m 2 q 2 ,q 2 ef ) and C Q (m 2 <? 2 ,<7 2 ef ) denote the DM scattering 
cross section and solar capture rate in the «; q 2 ,q 2 ef case, respectively. One may see from 
Fig. |4] (left panel) that <x„ will remarkably increase as increases when ~ q K f = 0.1 GeV. 
For < 0.01 GeV and > 0.2 GeV, the predicted cr n is insensitive to m^. These features can 
be easily understood from the forms of Fp M (g, v, m^). For the light DM mass m D = 10 GeV, our 
numerical results show that 6 kinds of dependent DM form factors can produce the similar 
curves. As shown in Fig. |4](right panel), the predicted C Q approaches to a constant as well as cr n if 

> 0.2 GeV. For < 0.2 GeV, C Q can usually decrease as increases. We find that the q 2 m^, 
q 4 m$ and V 2 q 2 m ( p cases have the minimums around « 0.04 GeV for C Q . In fact, the DM solar 
capture rates with a fixed cr n in the q 2 m,p, q 4 m^, and V 2 q 2 m l / ) cases are the monotone decreasing 
functions of m$. Therefore the minimums arise from the monotone increasing cr n . In terms of the 
results in Fig. |4] the DM solar capture rate in the m$ dependent scenario will quickly move from 
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the dashed line to the corresponding color solid line as increases in Fig. |3] When m$ > 0.2 
GeV, the dependent scenario will approach to the heavy mediator mass scenario. 



V. DISCUSSIONS AND CONCLUSIONS 

So far, we have used the usual Helm nuclear form factor for F 2 N (q) in Eqs. (OQ) and (fl4l) to calcu- 
late cr n and C©. In fact, the exact F 2 N {q) contains the standard SI nuclear form factor (Helm nuclear 
form factor) and an important correction from the angular-momentum of unpaired nucleons within 
the nucleus for the NR operators <9 3 and 4 in Eq. © IU8I1 . The correction is comparable with the 



standard SI nuclear form factor for nuclei with unpaired protons and neutrons when mo > m^. By 



use of the relevant formulas in Appendix A of Ref. H18H . we numerically calculate this correction 
contribution to the XENON 100 and XENON 10 experiments and find that it can be neglected for 
our analysis about cr n . In the previous sections, the predicted C© dominantly arises from the contri- 
butions of 'H, 4 He, 12 C, 14 N, 16 0, Ne, Mg, Si, S and Fe. Since these elements or dominant isotopes 
have not the unpaired protons and neutrons within the nucleus, our numerical results about C© are 
not significantly changed. 

In conclusion, we have investigated the SI momentum and velocity dependent DM-nucleus 
interactions and discussed their effects on <x„ and C©. In terms of the NR effective theory, we 
phenomenologically discuss 10 kinds of momentum and velocity dependent DM form factors 
Fp M (^r, v). Using these DM form factors, we have recalculated the corresponding upper limits on 
cr„ from the XENON100 and XENON10 experimental results. Each upper limit on cr„ can be 
used to calculate the corresponding maximal DM solar capture rate C©. Our numerical results 
have shown that the momentum and velocity dependent DM form factor cases can give larger 
DM annihilation signals than the usual contact interaction case for almost the whole parameter 
space. The light mediator mass scenario can usually produce the larger C© than the corresponding 
heavy mediator mass scenario except for the q 4 and 1 cases. On the other hand, we have also 
deduced the Super-Kamiokande's constraints on C©/2 for 8 typical DM annihilation channels 
with the equilibrium assumption Y A = C /2. In contrast to the usual contact interaction, the 
Super-Kamiokande and IceCube experiments can give more stringent limits on cr n than the latest 
XENON 100 experiment for several momentum and velocity dependent DM form factors when 
m D > 10 GeV. In addition, we have also considered 6 kinds of dependent DM form factors and 
analyzed their effects on <x„ and C©. We find that C© will quickly move from the dashed line to the 
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corresponding color solid line as increases in Fig. |3] When > 0.2 GeV, the dependent 
scenario will approach to the corresponding heavy mediator mass scenario. 
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